Constraints for weakly interacting light bosons 
from existence of massive neutron stars 
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Theories beyond the standard model include a number of new particles some of which might be 
light and weakly coupled to ordinary matter. Such particles affect the equation of state of nuclear 
matter and can shift admissible masses of neutron stars to higher values. The internal structure 
of neutron stars is modified provided the ratio between coupling strength and mass squared of a 
weakly interacting light boson is above g 2 /p 2 ~ 25 GeV~ 2 . We provide limits on the couplings 
with the strange sector, which cannot be achieved from laboratory experiments analysis. When the 
couplings to the first family of quarks is considered the limits imposed by the neutron stars are not 
more stringent than the existing laboratory ones. The observations on neutron stars give evidence 
that equation of state of the /3-equilibrated nuclear matter is stiffer than expected from many-body 
theory of nuclei and nuclear matter. A weakly interacting light vector boson coupled predominantly 
to the second family of the quarks can produce the required stiffening. 

PACS numbers: ll.10.Kk, 14.70. Pw, 26.60.Kp, 97.60.Jd 



Dark energy explains the accelerating expansion of the 
Universe. The density of dark energy po ~ 3.8 keV/cm 3 
may correspond to a fundamental scale \r> = Pjj^ ~ 
8.5x 10" 5 m [HEIEHS]. Theoretical schemes with ex- 
tra dimensions suggest modifications of gravity below A r> 
and a multitude of states with masses above 1/Arj very 
weakly coupled to members of multiplets of the standard 
model. Scales significantly below \r> represent the inter- 
est for supersymmetric extensions of the standard model 
which include generally a number of new particles, such 
as the leading dark matter candidate neutralino. Typi- 
cally, new particles are expected with masses above sev- 
eral hundred GeVs or even higher. However, light parti- 
cles may exist also, such as a neutral very weakly coupled 
spin-1 gauge [/-boson [5] that can provide annihilation of 
light dark matter and be responsible for the 511 keV line 
observed from the galactic bulge [SI [7] . 

Deviations from the inverse-square Newton's law are 
parametrized often in terms of the exchanges by hypo- 
thetical bosons also. Constraints on the deviations from 
Newton's gravity have been set experimentally in the 
sub- millimeter scale [5J [TOJ El DU CE3] and down to dis- 
tances ~ 10 fm where effects of light bosons of extensions 
of the standard model can be expected [TH[T51HHl[T71ll8j . 
Constraints on the coupling constants from unobserved 
missing energy decay modes of ordinary mesons are dis- 
cussed in Ref. [H)] . 

Bosons with small couplings escape detection in most 
laboratory experiments. However, bosons interacting 
with baryons modify the equation of state (EOS) of nu- 
clear matter. Their effect depends on the ratio between 
the coupling strength and the boson mass squared, so a 



weakly interacting light boson (WILB) may influence the 
structure of neutron stars even if its baryon couplings are 
very small. 

The effect of a vector boson on the energy density of 
nuclear matter can be evaluated by averaging the corre- 
sponding Yukawa potential: 

If g 2 e~v r 

El = 2 / dxi<2x 2 p(xi) — — — p(x 2 ), (1) 

where /o(xi) = p(xi) = p is the number density of ho- 
mogeneously distributed baryons, r = |x 2 — XjJ, g is the 
coupling constant with baryons, and p is the boson mass. 
A simple integration gives 
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where V is the normalization volume. 

The coherent contribution to the energy density of nu- 
clear matter from vector WILBs should be compared to 
that from the ordinary w-mesons. In one-boson exchange 
potential (OBEP) models, the nucleon-nucleon repulsive 
core at short distances r < b — 0.4 fm is attributed to u>- 
meson exchanges. Respectively, the w-meson plays a fun- 
damental role in nuclear matter EOS. In the mean-field 
approximation, the contribution of cj-meson exchanges to 
the energy has the form of Eq.Q, with g and p replaced 
by the ui- meson coupling g u and the mass p u . 

The NN interactions are described with g 2 //i 2 — 



175 GeV" 



The rclativistic mean field (RMF) model 



[2T] gives gli/pZi — 196 GeV . The compression modu- 
lus of nuclear matter K = 210 -7- 300 MeV is consistent 
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with gl/nl = 125 4- 180 GeV" 2 [12] . Stiff RMF mod- 
els use g 2 ,/ pt up to 300 GeV -2 [22]. If we wish to stay 
within current limits and do not want to modify the inter- 
nal structure of neutron stars qualitatively, as described 
by realistic models of nuclear matter, one has to require 
that vector WILBs fulfill constraint 

2 2 

% < % « 200 GeV" 2 . (3) 

A similar reasoning applies to scalar WILBs which 
have to compete with the standard cr-meson exchange. 
In OBEP models, the long-range attraction between nu- 
cleons is attributed to cr-meson exchanges. The contri- 
bution of the CT-mesons to the interaction energy has the 
form of Eq.Q, with g and p replaced by the cr-meson 
coupling g a and the mass p a . The sign of the contribu- 
tion must be negative because of the attraction. Also, 
p should be replaced by the scalar density. In RMF 
models, the cr-meson mean field decreases the nucleon 
mass. The effect depends on the ratio g 2 / p? also and 
produces an additional decrease of the energy at fixed 
volume and baryon number. The empirical values of the 
ratio g 2 /p 2 are 40 ~ 60% higher than those of the ui- 
meson (20j [2TJ [22l |23]. The internal structure of neutron 
stars is not modified significantly provided the coupling 
strength g and mass p of scalar WILBs fulfill constraint 

2 2 

^ < % « 300 GeV" 2 . (4) 

The deviations from the Newton's gravitational poten- 
tial are usually parametrized in the form 

\r( \ G-m 1 m 2 ( _ r/x \ . , 

V{r) = \l + a G e 'J. (5) 

The second Yukawa term can be attributed to new bosons 
with Gm 2 otG = ±5 2 /(47r) and A = 1/fi, where +/— 
stands for scalar/vector bosons and m is the proton mass. 

On Fig. 1 we show regions in the parameter spaces 
(g 2 ,p) and (ac,A) allowed for WILBs by the constraint 
([3]) . The constraint for scalar bosons is close t o C on- 
straints from other works [10] HU [El EH O [15] [16J EO 
[T8] are shown also. 

An increase of g (a decrease of p) of scalar WILBs in- 
creases the negative contribution to pressure, makes EOS 
of nuclear matter softer, makes neutron stars less stable 
against gravitational compression. The ratio g 2 / ' p? can- 
not be increased significantly above the limit Q, since 
the maximum mass of the neutron star sequence cannot 
be moved below masses of the observed pulsars. 

An increase of g (a decrease of p) of vector WILBs, 
conversely, increases the positive contribution to pres- 
sure, makes EOS of nuclear matter stiffer, makes neutron 
stars more stable against gravitational compression and 
drives the maximum mass of neutron stars up. 

In case of vector bosons, it is less obvious what kind of 
the observables confronts to high ratios g 2 / p 2 ■ 
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FIG. 1: (color online) Constraints on the coupling strength 
with nucleons g 2 /(47r) and the mass p (equivalently ac and 
A) of hypothetical weakly interacting light bosons: 1 are con- 
straints from Ref. (IS). 2 - from Ref. [TT], 3 - from Ref. [12] . 
4 - from Ref. [IS], 5 and 10 are constraints from low-energy 
n— 208 Pb scattering |16j and [14], respectively, 6 - from Ref. 
[17] . 7 - from Ref. [15], 8 and 9 are constraints from spec- 
troscopy of antiproton atoms [16], 11 and 12 are constraints 
from near-forward pn scattering for vector and scalar bosons, 
respectively [TS]. The axes are in the log 10 scale. The internal 
structure of neutron stars is not modified qualitatively pro- 
vided the boson coupling strengths with baryons and masses 
lie at g 2 /p 2 < 200 GeV" 2 beneath the highlighted area 13. 



Realistic models of nuclear matter are based on the 
nucleon-nucleon scattering data. They split into soft and 
stiff models according to the rate the pressure increases 
with the density. The soft models correspond to low max- 
imum masses of neutron stars ~ 1.6 M Q , while the stiff 
models give the upper limit around ~ 2.6 M . 

The problem on the softness of nuclear EOS has re- 
ceived new interest due the analysis of strange particle 
production in heavy-ion collisions. The data at different 
bombarding energies lead to the conclusion that EOS of 
nuclear matter must be soft at densities two to three 
times of the saturation density [2H US] [53] . Data on the 
transverse and elliptic flows in heavy-ion collisions sug- 
gest a soft EOS around the saturation, too [27] . 

Last years observations of pulsars with high masses 
have been reported. The most massive pulsars are PSR 
B1516+02B in the globular cluster M5 with the mass 
of 1.96+g;?| M and PSR J1748-2021B in the globu- 
lar cluster NGC 6440 with the mass of 2.74 ± 0.22 M 
|28j . The mass of rapidly rotating neutron star in the 
low mass X-ray binary 4U 1636-536 is estimated to be 
M = 2.0±0.1 M gS]. The mass and radius of the 
X-ray source EXO 0748-676 are constrained to M > 
2.10 ± 0.28 M and R > 13.8 ± 1.8 km [30]. The ob- 
servations on neutron stars suggest that EOS of the j3- 
equilibrated nuclear matter is stiff. 

The controversy between the conclusions on the soft- 
ness of nuclear matter as derived from the laboratory 
experiments and on the stiffness of the /3-equilibrated 
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nuclear matter as derived from the astrophysical obser- 
vations has been of interest since after the discovery of 
millisecond pulsars [3TJ[32] and earlier [3"3"] . 

Current models use to match EOS of neutron matter 
with a soft EOS at the saturation density and a stiff EOS 
at higher densities. Such models are in the qualitative 
agreement with laboratory and astrophysical data [34] . 

High densities provide favorable conditions for the oc- 
currence of exotic forms of nuclear matter: pion, kaon, 
and dibaryon condensates, quark matter. New degrees of 
freedom make EOS softer, pushing the maximum mass 
of neutron stars down. The recent astrophysical observa- 
tions seem to exclude the softest EOS e.g. based on the 
classical Reid soft core model [35] and make it problem- 
atic to accommodate the exotic forms of nuclear matter 
with masses and radii of the observed pulsars [3D] (see 
however [36]). 

The in- medium masses of vector mesons depend on the 
density. Assuming /x is a function of p and using Eq.Q, 
one may evaluate the w-meson contribution to pressure: 

A positive shift of the w-meson mass decreases the pres- 
sure and leads to a softer EOS, whereas a negative shift 
leads to a stiffer EOS. The data on the dilepton pro- 
duction in heavy-ion collisions do not give evidence for 
significant mass shift [37], so the observed stiffness of the 
/3-equilibrated nuclear matter can hardly be attributed 
to in-medium modifications of the vector mesons. 

The realistic models of neutron matter discussed in 
Ref. [34] neglect hyperon channels e.g. reactions £~ — > 
n + e+v e . In RMF models gHEHlES], the /^-equilibrium 
of hyperons drops the limiting mass by 0.5 ~ 0.8 M Q . 
This result is in accord with hypernuclear data and other 
recent calculations [101 EH H2] • The inclusion of the /3- 
cquilibrium for all baryons brings difficulties in reproduc- 
ing the observed masses of neutron stars. 

Coming back to vector WILBs, we see that their exis- 
tence is desirable to provide additional stiffening of the 
/3-equilibrated nuclear matter. 

The Compton wavelength of WILBs is assumed to be 
greater that the radius of nuclei e.g. l/^i > JJ w 7 fm 
ss (30 MeV)" 1 for the lead. The contribution of WILBs 
to the binding energy of nuclei then equals ~ A 2 g 2 /R 
like for photons. Since g 2 /(47r) is much smaller than the 
fine structure constant, the effect of WILBs on nuclei 
is negligible. Above ~ 10 2 MeV the coupling constant 
of WILBs is close to unity, so WILBs there are neither 
weekly interacting nor light. 

WILBs thus do not modify observables in laboratory 
experiments on hypernuclear physics, nuclear structure 
and heavy-ion collisions, since their baryon couplings are 
very small. The characteristic scale of the parameters of 
these particles is fixed by the upper limit 

The mass-radius relations for non-rotating neutron 
stars are shown on Fig. 2 for four values of the ratio 
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FIG. 2: (color online) Mass of non-rotating neutron stars 
as a function of radius: 1 - RMF model of hyperon mater 
with the compression modulus K — 300 MeV [55]; 2 - the 
same as 1 including a flavor-singlet vector WILB coupled to 
baryons with g 2 /p 2 = 25 GeV" 2 (1/8 of the limit (3)); 3 
- the same as 2 with g 2 / p 2 — 50 GeV -2 ; 4 - the same as 
2 with g 2 /p 2 = 100 GeV" 2 . The highlighted area within 
M = 1.96to ?2 M shows the mass constraint from PSR 
B1516+02B. The neutron star sequences should cross the ro- 
tation speed limit curves shown for pulsar PSR B1937+21 
with the rotation frequency of v = 642 Hz [51] and the neu- 
tron star XTE J1739-285 showing X-ray burst oscillations 
with frequency of v = 1122 Hz [52]. The mass-dependent 
lower bound on radii of neutron stars determined from the 
blackbody radiation of RX J1856. 5-3754 is shown. The dot- 
ted straight lines z = 0.1-^0.6 indicate the red shift at surfaces 
of neutron stars. The red shift of z = 0.35 measured for EXO 
0748-676 constrains the radii of neutron stars by R > 12 km 
and, respectively, masses [301 . 



g 2 /^ 2 = 0, 25, 50 and 100 GeV" 2 of a flavor-singlet vec- 
tor WILB. At densities below pdrip = 4.3 x 10 11 g/cm 3 
the matter represents an atomic lattice. WILBs do 
not modify properties of nuclei and the Baym-Pethick- 
Sutherland EOS [43], accordingly. At densities pdrip < 
P ;$ Pnuci — 2.8 x 10 14 g/cm 3 , atomic lattice coexists 
with neutron liquid. The matter at pdrip < p ^ Pnuci is 
described by the Baym-Bethe-Pethick EOS [44]. Above 
Pnuci 1 nuclei dissolve and the matter is described by the 
/3-equilibrated hyperon liquid with the compression mod- 
ulus K = 300 MeV [52]. WILBs contribute to the energy 
density and pressure above pdrip: as described by Eqs.([2]) 
and ([6| with dp/ dp — 0, through the spatially extended 
nucleon and hyperon liquid components of the neutron 
star matter. The vector WILBs give equal contributions 
to the chemical potentials of the octet baryons and do not 
violate the chemical /3-equilibrium. [53] The inclusion of 
such vector bosons does therefore not change composition 
of the neutron star matter. 

The highlighted area at the upper left corner of Fig. 
2 excludes within general relativity the radii of neutron 
stars below the Schwarzschild radius. The causal limit 
excludes the area R < 3GM [35]. The rotation speed 
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limit curves are constructed using the modified Keplerian 
rate u max ~ 1045 (M/M ) 1 /2(io km/i?) 3 / 2 Hz, which 
accounts for the deformation of rotating neutron stars 
and effects of general relativity [46 . 

It is seen from Fig. 2 that, despite we selected EOS 
with the high compression modulus, the neutron star se- 
quence with <? 2 //i 2 = contradicts to the mass mea- 
surement of PSR B1516+02B. It gives a very low mass 
of the neutron star from the blackbody radiation radius 
constraint also, which confronts with the lower limit of 
~ 0.85 Mq for masses of protoneutron stars [35] . 

The value of <? 2 //i 2 = 200 GeV -2 gives the maximum 
mass slightly above 3.0 M . However, the neutron star 
sequence does not cross the rotation speed limits, while 
the red shift remains always below z = 0.35. The upper 
bound ([3]) is thus critical for the internal structure of 
neutron stars. [54] 

The vector WILBs increase the minimum and max- 
imum mass limits and radii of neutron stars and are 
able to bring in the agreement models of hyperon mat- 
ter which are soft with the astrophysical observations 
on neutron stars which require a stiff EOS. The ratio 
g 2 / fji 2 ~ 50 GeV -2 might be reasonable. Such a value, 
however, clearly contradicts to the laboratory constraints 
shown on Fig. 1 in the entire mass range fj, = 10 -9 to 
10 2 MeV. 

The in- medium modification of masses of vector bosons 
modify EOS. Vector WILBs can be compared to the u>- 
meson where \6 fj, u \ / < 0.1 above the saturation den- 
sity [37] ■ A vector WILB mass shift can be estimated 
as S/i 2 ~ g 2 /g 2 2/J, u d fx u . The in-medium modification is 
small provided \S[i 2 \ < fi 2 i.e. g 2 7 'fi 2 < 10 3 GeV -2 , so in 
the region of interest Q holds for the vacuum masses. 

The laboratory constraints shown on Fig. 1 do not ap- 
ply to WILBs coupled to hyperons. A vector WILB cou- 
pled predominantly to the second family of the quarks 
makes hyperon matter EOS stiffer also. It contributes 
differently to chemical potentials of the octet baryons 
and suppresses the hyperon content of the neutron star 
matter due the additional repulsion. One can expect the 
ratio g 2 / ' (x 2 should be close to or higher than that esti- 
mated above (~ 50 GeV -2 ). In such a scenario, nuclear 
matter without hyperons can be treated as reasonable ap- 
proximation for the modeling structure of neutron stars 
in the /3-equilibrium also e.g. on line with Ref. [34J where 
models with the blocked hyperon channels are shown to 



be in the qualitative agreement with the laboratory and 
astrophysical constraints. 

Gauge bosons interact with the conserved currents 
only, but flavor is not conserved. A WILB coupled to 
the second family of the quarks cannot be a gauge boson, 
so it does not arise naturally in the current theoretical 
schemes. Here, we do not have a goal whatsoever to go 
beyond the phenomenological analysis. 

Hypernuclear data restrict NY potentials, whereas the 
interaction between hyperons YY is not known experi- 
mentally. The stiffness of the hyperon matter might also 
be attributed to the </>(1020)-meson exchange, whose cou- 
pling to the nonstrange baryons is suppressed according 
to the Okubo-Zweig-Iizuka rule (see, however, [50]). 

Summarizing, we have assumed the existence and de- 
rived constraints for a new boson that couples to nuclear 
matter. Such a particle contributes, by its coherent force 
among nuclear constituents, to a modified EOS and af- 
fects the structure of neutron stars. The neutron stars 
exclude scalar bosons with the coupling strengths and 
masses above the line 13 on Fig. 1, whereas in a nar- 
row band below it and above a vector boson coupled to 
quarks of the second family could modify the EOS in 
a direction favored by the observed masses and radii of 
neutron stars. The astrophysical constraints in the non- 
strange sector are less stringent than the most accurate 
laboratory ones. They are unique, however, for scalar 
WILBs in the strange sector. The region of validity of 
the astrophysical constraints extends from A ~ 10 fm 
to about 10 km. Detailed studies of manifestations of 
new bosons in astrophysics, physics of neutron stars, and 
hadron decays to energy missing channels can shed more 
light on the existence of WILBs and their possible effect 
on the structure of neutron stars. 
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